Even-odd effect in Andreev Transport through a Carbon Nanotube Quantum Dot 
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We have measured the current(/)-voltage(V) characteristics of a single-wall carbon nanotube 
quantum dot coupled to superconducting source and drain contacts in the intermediate coupling 
regime. Whereas the enhanced differential conductance dl/dV due to the Kondo resonance is 
observed in the normal state, this feature around zero bias voltage is absent in the superconducting 
state. Nonetheless, a pronounced even-odd effect appears at finite bias in the dl/dV sub-gap 
structure caused by Andreev reflection. The first-order Andreev peak appearing around V — A/e 
is markedly enhanced in gate-voltage regions, in which the charge state of the quantum dot is odd. 
This enhancement is explained by a 'hidden' Kondo resonance, pinned to one contact only. A 
comparison with a single-impurity Anderson model, which is solved numerically in a slave-boson 
meanfield ansatz, yields good agreement with the experiment. 

PACS numbers: 74.78.Na,74.45.+c,73.63.Kv,73.21.La,73.23.Hk,73.63.Fg 



There is a growing interest in the exploration of 
correlated charge transport through nanoscaled low- 
dimensional systems involvingboth superconductors and 
normal metals [jj, 12, y, Ij, la |fj| . The penetration of the 
pair amplitude A from a superconductor (S) into a nor- 
mal metal (N), the proximity effect, is a manifestation of 
correlated charge transport mediated by Andreev pro- 
cesses taking place at the S-N interface 0] and lead- 
ing in S-N-S junctions to the Josephson effect Q and 
sup-gap current peaks due to multiple Andreev reflec- 
tion (MAR) @. The superconducting proximity effect 
has been studied in great detail in the mesoscopic size 
regime of diffusive, but phase coherent conductors [Tol ] . 
Andreev transport has also been the key quantity in ex- 
periments elucidating charge transport in single atom 
contacts [ll[ . On the other hand, Andreev transport 
through a quantum dot coupled to superconductors, is 
just emerging now 0, [ID, [lj, 0- If the dot is weakly 
coupled to the leads, Andreev processes are suppressed 
by the charging energy U of the dot 0, O ■ If the dot 
is sufficiently small, size quantization takes place, form- 
ing a quantum dot (QD) with discrete eigenstates ('lev- 
els') at energies Esi}- Transport then occurs through 
individual levels Q. Since the level 'positions' Eu\, 
and sometimes also the coupling strengths of the lev- 
els to both source and drain contacts I\2, can be tuned 
through gate voltages, a physically tunable model system 
of the Anderson 'impurity problem' is realized. With 
one electron on the QD (half- filling), a many-electron 
ground-state forms, involving both the dot-state and con- 
duction electrons from the leads in an energy window 
given by the the Kondo temperature Tk 0, 0]- In 
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this Kondo regime, which can be observed if ^1,2 is not 
too small, a resonance pinned at the Fermi energy of 
the leads forms (Kondo resonance). If superconduct- 
ing contacts are used instead of normal ones, the ad- 
ditional pair-correlation in the leads competes with the 
Kondo correlations on the QD 0, @, 0, EE E3, 0. 
It has been found recently in experiments using carbon 
nanotubes (CNTs) as QDs, that there is an interesting 
cross-over occurring at ksTx * A. If A > fc^T^, the 
Kondo correlations are suppressed, whereas they persist 
in the opposite regime, opening a highly conducting chan- 
nel for the Josephson effect [2|, |l2j. CNTs are ideally 
suited for the realization of such systems, because CNTs 
can a) act as well controlled QDs in different transport 
regimes [24[, including the Kondo regime [251 ]. and b), dif- 
ferent types of contacts can be realized, including super- 
conducting ones 0, i, [H, [H H 0. Similar physics 
can be addressed with QDs fabricated in semiconducting 
nanowires contacted to S contacts [26| . 

We here report on measurements of the non- 
equilibrium (finite-bias) transport through a single-wall 
carbon nanotube QD with S contacts in the most in- 
teresting regime of intermediate coupling, where Kondo 
correlations are of similar magnitude as superconducting 
ones. We have found a pronounced even-odd effect in the 
MAR structure. 

Single walled carbon nanotubes were grown by chem- 
ical vapor deposition on highly-doped Si wafers with a 
400 nm layer of thermal oxide on top, using Fe parti- 
cles as catalysts [27|, l28[. Individual single- wall carbon 
nanotubes (SWCNTs) were localized with a scanning 
electron microscope and contacted to supercondcuting 
source and a drain electrodes using e-beam lithog- 
raphy, see Fig. 1. The evaporated contacts consist 
of a Ti(5nm)/Al(100nm)/Ti(10nm) tri-layer, where Ti 
serves as an adhesion and cap-layer. Al is the ac- 
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FIG. 1: (a) Illustration of a SWCNT contacted by supercon- 
ducting (Al) source (1) and drain (2) electrodes. fj,i t 2 are 
respectively the electrochemical potentials, Nip the density- 
of-states, I\2 the life-time broadenings due to the coupling of 
the QD state to the electrodes, A the superconducting gap pa- 
rameter, V the applied source-drain voltage and V g the gate 
voltage, (b) Device geometry, showing the evaporated tri- 
layer, consisting of a Ti adhesion layer, the superconducting 
Al film and a Ti cap-layer, (c) shows an actual device and 
(d) illustrates possible processes that lead to a subgap cur- 
rent. Shown is in solid a first order Andreev process and 
dashed a second order one. In the first (second), two (three) 
quasiparticle (electrons e and holes h) are involved. 



tual superconductor with a bulk critical temperature of 
T c = 1.2 K. In its thin film form, we rather measure a 
T c of 0.9 K, which corresponds to a BCS gap-parameter 
A = 1.76k B T c of 0.135 meV. We drive the Al contacts 
into the normal state by applying a small perpendicular 
magnetic field of B = 0.1 T. The substrate is contacted 
to a third terminal in order to establish a backgate. We 
measure the differential source(l)-drain(2) conductance 
G := dl/dV as a function of source-drain V and gate- 
voltage V g . This is achieved by superposing an ac- voltage 
V ac = 10 /iV on V and measuring the corresponding ac 
current with a current-voltage converter and a lock-in 
amplifier. Several devices were fabricated and tested at 
room temperature and at 4.2 K. Here, we focus on a par- 
ticular interesting device which we selected and measured 
in a dilution refrigerator. This device has been studied 
over a large V g window and displays single-electron charg- 
ing with addition energies in the range of 2 ... 5 meV. In 
the following we will focus on a confined gate-voltage 
regime. 

Fig. 2 shows the main measurements we will be focus- 
ing on in the following. Visible is in (a) a dl/dV plot 
in the normal state (n-state) and the corresponding one 
in the superconducting state (s-state) in (b). In the re- 
state a sequence of larger and smaller Coulomb blockade 
(CB) diamonds are seen (dashed lines), corresponding to 
a sequence of nearly equidistantly spaced levels on the 
SWNT-QD, which are filled sequentially. The number of 
electrons on the dot therefore alternates in the ground- 
state between odd and even [29j. It is also seen that the 
conductance G — dl/dV around zero bias is suppressed 




FIG. 2: Differential conductance dI/dV(V, Vg) plot as a 
function of bias V and gate voltage V g of a SWCNT-QD 
with superconducting contacts in the normal (a) and super- 
conducting state (b). In (c) we show the linear conductance 
G(T) measured as a function of temperature T in the middle 
of charge state 3. The curves in (d) correspond to the one 
overlaid on the dl /dV(V, V g ) plots. 

and featureless in the even valleys, but is increased as- 
suming structure in the odd ones. In the CB diamond 
labelled 3, there is a pronounced peak at V = 0, sug- 
gesting the appearance of a Kondo resonance. Indeed, 
the dependence of the linear G(T) on temperature T 
(Fig. 2c), measured in the middle of this valley, follows 
the expected dependence (30j with a Kondo temperature 
of Tk — 0.75 K. In the other odd valleys, the Kondo res- 
onances are split by re 0.1 meV 31]. The origin of this 
splitting is at present not known, but could be due to ex- 
change with ferromagnetic catalyst particles or another 
tube (if the present one is a small bundle or contains 
more than one shell) [32| . 

The n-state data can be used to deduce a number of 
parameters. The source, drain and gate capacitances 
are G\^ ~ 50, 100 aF and C g ~ 4aF, leading to a gate- 
coupling a = Cg/Cs of ~ 0.026, where Cs = C\ + 
C*2 + Cg. The charging energy U = e 2 /Cs and the 
level spacing 5E are in the range of 0.7. . . lmeV and 
1.4 . . . 1.8 meV, respectively. Whereas this SWNT-QD is 
nearly symmetric in its electrostatic coupling, it is quite 
asymmetric in its electronic one. The total level broad- 
ening amounts to T = T± + T2 re 0.2 meV with an asym- 
metry of ri/r2 ~ 50. This asymmetry is deduced from 
the measured current peaks in dl/dV at the border of 
the CB diamonds at finite bias and is in agreement with 
the reduced low temperature zero-bias G(0) of the Kondo 
ridge 3, amounting to G(0) ~ 0.1e 2 /h. 

Looking next at the s-state, we see that the major 
changes in the dl/dV are confined to a voltage band 
of — 0.26meV< V < 0.26 meV, corresponding to ±2A. 
Above 2A, i.e. |V| > 2A/e, quasiparticle current is pos- 
sible and the main modification is caused by the peak 
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FIG. 3: (a) measured dl /dV versus temperature T in the even 
valley 4 of Fig. 2a. The 2 A and A peak positions E pBa k (T) are 
shown in (b) together with a BCS T-dependence of A. In (c) 
we show the result of a model calculation based on tunneling 
between two superconductors. Good agreement is obtained if 
the BCS-DOS is broadened by an amount corresponding to 
0.2Ao. The illustration in (d) is our proposal to understand 
the appearance of the strong A feature in the odd valleys. 
A Kondo resonance persists on one electrode side only (here, 
drawn as a thin curve on the left side), pinned to that chem- 
ical potential. Due to the gaped DOS, this resonance splits 
leading to an enhancement of the A feature. 

in the superconducting density-of-state (DOS) [!, [TlT |. 
leading to a peak-like feature in dl/dV. Due to the 
gap in the DOS, first order processes are however sup- 
pressed below 2A. Depending on the transmission prob- 
ability [U 0, E3] , higher order Andreev processes can 
contribute. The first order Andreev process, for example, 
results in a peak- like structure in the vicinity of A. Due 
to the higher order, the first Andreev peak and all sub- 
sequent ones are usually smaller than the quasiparticle 
one. Both the dominant 2A and the smaller A peak, as 
well as the suppressed G in the s-state are best visible in 
the middle of an even charge state (even valley), see e.g. 
curve labelled 4 in Fig. 2d. In contrast, in the odd charge 
states, the 2A feature is not present or does not appear 
at 2A. Starting to view the data from large bias voltage, 
the first peak appears closer to A rather than 2A, with a 
preceding negative dl/dV (NDR), see curve labelled 3 in 
Fig. 2d. Hence, there is a striking even-odd asymmetry 
in the finite-bias dl/dV features in the s-state which is 
not caused by the CB resonance at the edges of adjacent 
charge states, where the situation is expected to be more 
complex [l(| E3 • The even-odd alteration of the MAR 
structure suggest a relation to Kondo physics. To model 
this, we first extract important parameters from an anal- 
ysis of the data in the middle of an even valley where 
Kondo correlations are absent. 

Fig. 3a-c discusses the temperature dependence of 
dl/dV in the middle of the even charge state; (a) shows 
the measurement taken in valley 4 of Fig. 2 and (c) 



is a model calculation. The experiment displays pro- 
nounced quasiparticle current peaks at E pea k = ±2A, 
and weaker MAR peaks at ±A. The evolution of 
E^f ak (T) with temperature T are shown in Fig. 3b to- 
gether with an approximate BCS gap function E^ ak = 

Kj A A tanh(1.74 A /T c /T - 1), where we used the BCS 
value for Ao = 1.76fcsT c , which amounts to 0.135meV 
for a T c of 0.9 K. We then obtain K 2A = 2.0 and 
A'a = 1.15 for the two peaks. The slight increase of Ka 
above the expected value of 1 is not unusual. We will be 
using the value A = 0.13 meV as the zero-T gap param- 
eter in the following. The relevant parameters expressed 
in units of A are then: U = 5 ... 8, SE = 10 ... 14, 
r « 1.5 and Tk ~ 0.5. The zero-bias peak in Fig. 3a, ap- 
pearing at intermediate temperatures, can be explained 
by direct tunneling of quasiparticles thermally activated 
across the gap. 

The good agreement with the BCS relation of the peak- 
positions motivates the modelling of the dl/dV using the 
BCS-DOS in the leads. Although we would have to use a 
theory describing resonant tunneling between two super- 
conductors, such as the one from Levy Yayati et al. [l6|, 
a simple tunneling picture suffices [33|], because the res- 
onant levels are far away from the electrochemical po- 
tential of the electrodes in the middle of a charge state. 
The subgap current is treated in the same approximation 
using [ll|. To obtain a reasonable fit, the BCS-DOS has 
been convoluted with a Gaussian of width ij. A quite 
good agreeement is found with r] = 0.2 A . 

We now turn our attention to the odd charge states. 
We point out, that the zero-bias high-G Kondo 'ridge', 
which is associated with the Kondo resonance and visible 
in the n-state, is not seen in the s-state. This is consistent 
with a Kondo temperature Tk that is smaller than A, i.e. 
Tk = 0.5 • A @, [lj]. Although the Kondo resonance is 
not visible in the s-state, we suspect it to be responsible 
for the even-odd asymmetry of the A-feature in the s- 
state. 

In the Kondo regime, the single spin on the QD in 
the odd state is screened by exchange with conductance 
electrons from the leads. If the quasiparticles are bound 
in Cooper pairs in the leads, the Kondo temperature is 
renormalized, assuming a smaller value T K . This renor- 
malization is sensitive on the actual parameters Ti^ and 
A. Due to the asymmetry, it may happen that a Kondo 
resonance with a reduced width forms on the contact with 
the larger T, whereas on the other one it is suppressed. 
This is illustrated in Fig. 3d. 

We modelled this scenario considering a single-level 
Anderson Hamiltonian with interaction U = 5 . . . 10 Ao 
coupled to source and drain contacts. The experimen- 
tally deduced T's, including the strong asymmetry were 
used. The calculation is based on a slave-boson mean- 
field treatment of this interacting problem [2l[ . The re- 
sult of the comparison is shown in Fig. 4: (a) corresponds 
to the n-state and (b) and (c) to the s-state. Despite this 
simple model, the agreement is surprisingly good. It is 
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FIG. 4: Comparison of the dl/dV in the odd valley labelled 3 
in Fig. 2a with a model caluclation based a single-level Ander- 
son model with interaction U, that is evaluated by a mean- 
field slave-boson ansatz. The n-state is shown in (a), the 
s-state in (b) and (c). e denotes the level position. The ex- 
periment corresponds to U — 7 — 8 (in units of Ao), whereas 
the calulation was done for U = 5 and U = 10, where the 
latter is shown here. The upper solid (dashed) curve in (d) is 
cross-sections in the s-state at e = 0, taken from the experi- 
ment (theory). 

remarkably good in the normal state, shown in Fig. 4a. 
In the s-state, the dominance of the A-like feature in the 
odd valley is clearly present, as is a similar cross-over 
from odd to even filling. There are also some differences: 



in the experiment the A-feature bends to larger ^-values 
in the middle of the odd state, whereas this feature is 
rather flat in the calculation. 

In conclusion we have discovered a pronounce even- 
odd effect in the (multiple) Andreev structure in trans- 
port through a QD with superconducting contact. This 
effect is driven by a Kondo resonance pinned to one con- 
tact only and defines a new regime. Whereas a high 
conductance channel from source to drain, driven by 
Kondo correlations, persists in the supercondcuting state 
if Tk 3> A, this channel is greatly suppressed in the op- 
posite limit. In the intermediate regime Tk ~ A, and 
in particular for asymmetric dot-electrode couplings, the 
(partial) Kondo-scrccning of the 'impurity' spin may oc- 
cur on one electrode only. It would be interesting to 
explore the 'robustness' of this feature in model calcula- 
tion and to fabricate similar QDs with tunable electrode 
couplings. 
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